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Abstract 

We have carried out an ab-initio study of a-Sn (111), with the aim of pre- 
dicting and understanding its structure, reconstructions, and electronic states. 
We consider a variety of structural possibilities, and optimize them by moving 
atoms according to Hellmann-Feynman forces. Our results indicate that the 
unreconstructed surface is highly unstable, while a variety of reconstructions 
compete for the true ground state. Extrapolated trends from diamond to Si to 
Ge are well borne out, with a 2 x 1 vr-bonded chain reconstruction prevailing 
in the absence of adatoms, and a c(4 x 2) or (2 x 2) basic adatom-restatom 
unit reconstruction otherwise. Accompanying surface bucklings are in both 
cases larger than in Si and Ge, with consequently large ionic charge transfers 
predicted. Search for a /?-Sn-like metallic state of the surface turned out to 
be inconclusive. 
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I. INTRODUCTION 



Essentially all semiconductor surfaces are reconstructed. As is well known, this can be 
traced back to the strong covalency of their bulk band structure. In fact, the uniqueness of 
covalent systems, as opposed to, for example, regular metals, is that bulk covalent bonding 
is very directional and relatively inflexible. Consequently, it generally costs too much energy 
to locally rearrange the electron wave functions at the surface (as required by the presence 
of broken bonds), without some major lattice readjustments. 

In the group-IV insulators and semiconductors, known to crystallize in the diamond 
structure, the best studied face by far is the (111), where well documented results are 
available for diamond. Si and Ge. Both experiment and theory indicate that a variety of 
different reconstruction mechanisms are realized. A summary is presented in Table |. One 
can see that all (111) surfaces possess a vr-bonded chain 2x1 reconstructed state. As a 
true ground state, however, vr-bonded chains only prevail in diamond, possibly reflecting 
the difficulty of obtaining well annealed adatoms on this surface. On Si and Ge (111), the 
adatom reconstructions, (7 x 7) [Ref. |l| and c(2 x 8) [Ref. 0] respectively, have lower energy, 
whereas the 2x1 vr-bonded chain reconstruction remains well-defined, but energetically 
metastable. 

In all cases: diamond. Si and Ge (111), the metallic character of the ideal, symmet- 
rical vr-bonded chain reconstruction! is removed by some accompanying distortion. This 
distortion can consist either of a slight dimerization of the chain, as is probably the case in 
diamondi, or of a vertical buckling, as in SiiH3 and gM. Although experimental buckling 
data for diamond and Ge (111) are not yet available, ab-initio calculations exist for these 
geometrical distortions. The size of the predicted chain bucklings exhibits an interesting 
trend, namely it is zero in diamondi, about 0.4 A in siliconi (where experimental datai 
confirm this magnitude), and 0.8 A in germanium!. The absence of buckling in carbon has 
been ascribed to the larger intra-atomic Coulomb repulsion, which prevents the ionic charge 
transfer implied by the buckling. The decrease of Coulomb repulsion with increasing atomic 
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number is consistent with the trend towards a larger buckhng. 

When adatoms are present, as in annealed Si and Ge (111), then an adatom-rest 
atom reconstruction prevails. Interestingly, the adatom reconstruction again implies sur- 
face ionicityl^, with an electron charge transfer between the adatom and the restatom. 

Besides vr-bonded chains and adatom-restatom reconstructed states, a third kind of state, 
namely the displacive 2x1 buckled reconstruction, originally proposed by Haneman0, 
may also play a role. Although never found to be stable experimentally, and never stable 
theoretically either, it does make a brief appearance as a metastable transition state in the 
transformation path from ideal 1 x 1 to 2 x 1 vr-bonded chain, in Sii and Ge (lll)i, as found 
by ab-initio molecular dynamics. This does not seem to occur on C (111)@, probably again 
because the Haneman state is ionic, and thus very costly in carbon. 

Finally, the last entry in Table | considers, as a fourth type of reconstruction, the possible 
transformation of a thin surface layer from semiconducting to metallic. This possibility is 
suggested by the bulk phase diagrams of Si, Ge and Sn, where fully metallic phases com- 
pletely surround the semiconducting phases, under any combination of either high pressure 
or temperature. There is considerable evidence that surface metallization does in fact play 
a relevant role at high temperatureJll. The finite temperature behavior is however outside 
the scope of the present work. Although this transformation (as predicted, for example, in 
Ga0) does not appear to prevail at T=0 K in Si and Ge, it could still take place in a-Sn, 
which is closer to being metallic. 

Throughout the literature, only the (100) surface of a-Sn appears to have been consid- 
ered so fai§. This paper is devoted to what is therefore a first study of the a-Sn (111) 
surface. In particular, since there are no experimental data available, ab-initio calculations 
are required to compare the various possible reconstructions, and understand the trends 
relative to the other semiconductors just described. By extrapolating the trends of Table | 
we are led to expect the following: 

i) the 2x1 TT-bonded chain reconstruction should be present also in a-Sn(lll), with an even 
more pronounced buckling than in Ge; 



ii) This reconstruction will be in energetic competition with some adatom/restatom recon- 
structions, the latter being likely to prevail; 

iii) A surface metallization mechanism is also possible because of the presence of the bulk 
metaUic /3-Sn phase, energetically very close to the T — bulk ground state a-Sn. 

In this work we shall find that (apart from metallization, which turns out to be a harder 
problem) these trends are generally well borne out by ab-initio calculations. An adatom- 
restatom reconstruction is predicted to be the annealed ground state of ct-Sn (111). Al- 
though the exact periodicity cannot be predicted, our geometry and energetics based on an 
elementary c(4 x 2) or 2 x 2 unit should be fairly reliable. 

This paper is organized as follows. First of all, we calculate electronic structure, total 
energy and optimal geometry of bulk Sn, both in the a semimetallic phase, and in the 
P metallic phase. This is described in Section II. Next, in section III we use the same 
ingredients for a surface calculation in the usual slab geometry. We deal first with the 
ideal, unreconstructed a-Sn (111) surface, whose properties we study without and with 
surface relaxation. We find partly filled dangling bond states in the ideal surface electronic 
structure, which as usual suggest that this is not a stable surface. This is confirmed in 
Section IV, where a 2 x 1 displacive reconstruction is shown to lower surface energy without 
energy barriers. The more interesting 7r-bonded chain reconstruction is considered in Section 
V. It is found that a very strongly buckled 2x1 7r-bonded chain state exists for a-Sn (111), 
its energy being better than that of the Haneman state. Section VI is devoted to the 
alternative possibility of adatom/restatom reconstructions. We find this latter possibility to 
be energeticall ythe best, particularly in the c(4 x 2) local geometry. Lastly, in Section VII we 
briefiy describe an attempt at finding a stable metallic reconstructed surface ground state, 
which however turns out to be inconclusive. Section VIII summarizes the main conclusions 
of this work on a-Sn (111), confirming the trends extrapolated from diamond. Si and Ge, 
and stressing our predictions for future experimental work. 
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II. BULK PROPERTIES 

The zero temperature structure of bulk Sn is the a phase with a diamond latticelil of 
spacmg 6.483 A (at 90 K)lll. As is well known, if the temperature is raised to 286 K at 
atmospheric pressure, semimetallic a-Sn transforms into the fully metallic /5-Sn phase0. 
Under pressure, this transformation occurs at lower temperatures, extrapolating to T=0 K 
at a pressure of 5 Kbai0'0. Total energy LDA (local density approximation) calculationslH 
confirm that at zero temperature and zero pressure the (3 phase is energetically very close to 
the a phase, their energy difference being 44 meV/atom, with a predicted transition pressure 
to the /? phase of 8 Kbar. The relative success of LDA calculations in predicting very closely 
these extremely delicate differences provides confidence in extending this approach towards 
exploring surface properties. 

In the present ab-initio calculations for a-Sn (bulk and surface) and /3-Sn (bulk) we model 
the electron-ion interaction with a norm-conserving pseudopotentiallll in the Kleinman- 
Bylander formlll with s and p nonlocality, and use the standard LDA for the electron-electron 
interaction. We expand the Kohn-Sham orbitals in plane waves with an energy cutoff of 12 
Ry in all calculations (bulk and surface) for consistency. For the a-Sn bulk calculation, 10 
special k-points0 are sampled in the Irreducible Brillouin Zone (IBZ) of the diamond phase. 
For the metallic (3 phase, instead, we find that 160 k-points are needed to describe accurately 
the IBZ summations, with an additional gaussian spreading of 0.14 eVi. We find a lattice 
constant a = 6.446 A (exp. a = 6.483 aM) for a-Sn and a = 5.768 A (exp. a = 5.812 
Ai) for /?-Sn (for simplicity, we fixed the c/a ratio to the experimental value of 0.546). We 
also find that the /5-phase is disfavored with respect to the a-phase by 17 meV/atom, with 
an a /9 transition pressure of 5 Kbar in excellent agreement with experiment. We have 
tested that increasing the energy cutoff to 35 Ry and the k-point integration to 200 points 
in the IBZ of the /5-phase does not alter significantly these results. On the contrary, if the 
energy cutoff is less than 12 Ry, e.g. 10 Ry, the transition pressure becomes 2.5 Kbar and 
the energy difference drops to 9 meV/atom. 
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The electronic band structure of a-Sn is reported in Fig. ^ and that of /3-SnEl in Fig. 
^. As can be seen, in particular, the semimetallic character of the a-phase (zero gap at 
r) and the fully metallic character of the /5-phase are of course correctly borne out in our 
calculations. 

III. THE RELAXED, UNRECONSTRUCTED a-SN (111) SURFACE 

We study the surface properties of a-Sn (111) by a standard slab calculation. One of the 
two surfaces of the slab is frozen in its ideal geometry, together with the first three adjacent 
layers (a total of two rigid bilayers). The atoms belonging to the remaining layers are 
allowed to fully relax guided by the corresponding Hellmann-Feynman forces. Convergence 
is assumed when forces are less than 4 meV/A . The number of atomic layers in the slab 
is fixed to a total of 12 layers (not including the adatom layer when present). Each layer 
consists of either one, two, or four Sn atoms, in correspondence to choosing 1 x 1, 2 x 1, 
and 2 X 2 [or c(4 x 2)] surface cells, respectively. The number of "vacuum layers" is fixed 
to 6 (vacuum thickness ~ 11 A ), and the number of k points in the Irreducible Surface BZ 
(ISBZ) is chosen according to the size of the surface cell and its geometry, as described case 
by case. The initial atomic positions are chosen according to the calculated equilibrium bulk 
lattice spacing (oo = 6.446 A ). This guarantees that the forces on the atoms of the central 
layers are smaller than 4 meV/A for all the surfaces studied. In all cases, surface energies 
Esurf are given as Esurf = Esiab - NEb^ik - ^fur}^ where E^iah is the total energy of the slab, 
Ehuik is the energy per atom of bulk a-Sn as computed in section II [_Efe„/fc=-96.753 eV/atom] 
above, is the total number of atoms in the slab. Furthermore -E^ur/ = l^lfab'' ~ NEhuiklf^ 
is the energy of the frozen ideal surface. Here Elf^j}'' is the total energy of a slab with both 
surfaces rigid and ideal, computed using the same supercell geometry and k-point set as Egiab- 
We thereby have to repeat ideal surface calculation several times in different supercells and 
k-point sets corresponding to different surface reconstruction slab calculations in order to 
compare different reconstructed surface energies. The convergence of our results with respect 
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to the k-points sampling has been tested for each surface by increasing the k-point number 
after the atomic relaxation. Accordingly, we can estimate our overall energy resolution to be 
about 10 meV/(l x l)cell. We also constructed the electronic band structure of the relaxed 
surfaces, by using the Kohn-Sham eigenvalues of a 10-layer slab, obtained by removing the 
seven bottom layers of the simulation slab and replacing them by the inverted image of the 
topmost five relaxed layers. The reason for this procedure is threefold. First, we get rid 
of the undesired states related to the bottom rigid ideal surface. Second, we increase the 
effective symmetry of the supercell, thus decreasing the computational effort of the band 
structure calculation. Third, although the interaction of identical surface states belonging 
to opposite surfaces generally lifts their degeneracy (again an undesired effect), their average 
energy still corresponds, to first order, to the noninteracting value in the ideal case of an 
infinite slab. 

As a first case we have considered the (1 x 1) unreconstructed surface. Six special k- 
points in the hexagonal ISBZ have been used0. The electronic band structure of the ideal 
(1x1) surface is reported in Fig. ^. There are various surface states lying in the projected 
gaps. The surface states crossing the Fermi level inside the fundamental gap are clearly 
related to the presence of unsaturated surface dangling bonds. This feature is common to 
all the other group-IV insulator and semiconductors, and is responsible for the high grade 
of instability of the ideal (111) surface, as mentioned. The electron density corresponding 
to the surface state is analysed in Fig. ^, and reveals a high degree of surface localization 
as well clear dangling-bond character. 

As the next step, we allowed the surface, i.e. all atoms in the eight topmost surface layers, 
to relax, according to Hellmann-Feynman forces, so as to reduce surface energy. Despite the 
presence of unsaturated dangling bonds, the result of this energy minimization shows that 
surface atoms do not relax significantly (see Tabel y), with a surface energy gain of only 7 
meV/(l X l)cell, and a downward relaxation of the top layer of 0.02 A . We also checked 
that increasing the k-point number from six to eighteen changed none of the above results. 
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IV. HANEMAN (2 x 1) "BUCKLED ATOM" RECONSTRUCTION 

In the previous calculation, reconstruction was forbidden by symmetry. If the symmetry 
constraints imposed by the choice of a (1 x 1) cell are relaxed, the ideal surface is provided 
with a simple mechanism for the partial saturation of the dangling bonds. One such mech- 
anism, first proposed by Haneman0, consists of a simple in-out buckling of the topmost 
layer, resulting in a (2 x 1) displacive reconstruction. The inward motion of one surface 
atom implies an sp^ — sp'^ rehybridization, and a more p^-like dangling bond for that atom. 
The outward motion of the other atom, by contrast, causes dehybridization, and a more 
s-like dangling bond. Since in the atom Eg « Ep^, electrons will flow from the inward to 
the outward atom. This charge transfer empties and saturates respectively the two dangling 
bonds and therefore stabilizes the surface. This (2x1) reconstruction, although actually 
never observed, has been recently suggested to play the role of a transition state, a kind of 
"stepping stone" in the dynamical process leading from an unreconstructed state towards, 
e.g., a (2 X 1) TT-bonded chain statei. In order to study the possible occurence of such a 
buckled atom reconstruction on a-Sn (111), we have repeated our calculation in a larger 
(2 X 1) surface cell, using 4 special k-points in the rectangular (2 x 1) ISBzi. We find, 
indeed, that the surface spontaneously buckles as in the Haneman distortion, against which 
the ideal surface is therefore unstable, very much as Ge (111) doesi. The energy gain from 
the ideal to the optimal buckled geometry (Table ||), is about 0.22 eV/ (1 x l)ceU, measured 
relative to the energy of the ideal surface computed using the same (2x1) supercell and 
k-point set. The buckling of the top layer is found to be enormous, namely 1.23 A . The 
final atomic coordinates of the relaxed top five atomic layers are given in Table |IT|. 

V. (2 X 1) vr-BONDED RECONSTRUCTION 

The instability of the ideal (111) surface against (2 x 1) buckling is interesting, but prob- 
ably academic, except possibly in dynamics. In our pursuit of the true (111) ground state 
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surface structure, we consider next the surface reconstruction geometries which are experi- 
mentally observed in the other group-IV insulator and semiconductors. Here we consider the 
(111) 2x1 TT-bonded chain reconstruction. In the lack of any experimental and theoretical 
data, we arranged the surface atoms in the structure proposed by Pandeyli. Pandey's struc- 
ture can be obtained by simultaneously depressing one surface atom into the second layer 
and correspondingly raising one second layer atom, so that i) the number of dangling bonds 
per surface atom is the same as before; ii) all bond lengths are set to the bulk value. The 
atomic coordinates in Pandey's structure are given in Table Subsequently, we allowed 
the atomic positions to relax according to their ab-initio forces, using the same supercell 
and same k-point sampling as in the previous section. At convergence, the energy gain of 
this 2x1 reconstruction, with respect to the ideal surface, is 0.24 eV/ (1 x l)cell (Table 
HD, i.e., larger than that of the previous section. The optimal atomic coordinates of the top 
five atomic layers are given in Table The bond lengths within the vr-bonded chain are 
2.80 A , only 0.4 % longer with respect to the ideal bulk bond length (2.79 A ). The chain 
buckling (1.15 A ) is however very large if compared with that of Si or Ge (0.4 Ail, and 
0.8 AI, respectively), fully confirming the trends discussed in Section I and listed in Table 
|. The magnitude of this buckling is so large, that the whole 7r-bonded chain is now lying 
onto an essentially vertical plane (see Fig. ^). The large chain buckling is accompanied by a 
large electron transfer from the lowered to the raised chain atom, which largely saturates the 
raised atom dangling bond. This is clearly seen in Fig. ^, where the charge distributions of 
the highest occupied (panel (a)) and lowest unoccupied states (panel (b)) are shown. This 
charge rearrangement is accompanied by the opening of a large gap in the dangling bond 
surface state, as seen in the electronic structure of Fig. |[ Here the highest occupied (lowest 
unoccupied) surface states correspond to raised (lowered) atom dangling bonds respectively. 
Finally, since the alternative possibility of a dimerization of the chain was excluded in our cell 
due to symmetry constraints, we tried slightly dimerizing the initial, unbuckled 7r-bonded 
chain, and also the final, fully buckled one (thereby doubling the number of k-points in the 
ISBZ). However, we found these configurations to be energetically disfavored with respect 
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to the undimerized chain. 



VI. ADATOM / RESTATOM RECONSTRUCTIONS 

As the next hkely candidate for the reconstruction of a-Sn (111), we now consider the 
adatom / restatom reconstruction, known to be the stable mechanism for both Si and Ge. 
Although this reconstruction shows up with rather more complex surface unit cells, such 
as the (7 x 7) DAS model in Sii and the c(2 x 8) in Gel, the building block of this class 
of reconstructions is simple. It is based on the presence of one adatom every four (1 x 1) 
first-layer atoms. The adatom sits in a threefold site, saturating three first-layer atoms, 
leaving one (the restatom) unsaturated. Of the two available threefold sites, namely T4 (on 
top of a second- layer atom), and (hollow site), the adatoms prefer, at least in Si and 
GeHU, the T4 site. Moreover, the T4 site adatoms may still be arranged in a (2 x 2) or 
c(4 X 2) geometry. For instance in Ge (111) c(2 x 8), they are stacked in alternated (2 x 2) 
and c(4 x 2) cells. In this work we assume the pure (2 x 2) [or pure c(4 x 2)] T4 structure 
as the prototype adatom / restatom reconstruction, restricting for simplicity our analysis 
to this case only. As for the vr-bonded reconstruction case, we have no data to guide us. 
We start with an ideal T4 position for the adatom (see Table 0) such that the lengths of 
the adatom bonds with the three first-layer atoms are equal to the bulk bond-length. The 
relative positions of atoms are initially chosen according to the (2 x 2) [or c(4 x 2)] case 
of Table ^ [or Tabel ^|]. The fc-point sampling of this larger supercell is restricted to a 
single special point (mean-value point) of the hexagonal (2 x 2) ISBZE. We then calculate 
Hellmann-Feynman forces, and let the atoms relax to the equilibrium positions. We also 
find that increasing the number of special k-points of the relaxed surface to three and six 
does not change the surface energy within 5 meV/ (1 x l)cell. As before we must also repeat 
the rigid ideal slab calculation using this larger supercell and k-point sampling. As shown 
in Table the energy of the relaxed adatom-reconstructed surface turns out to be lower 
than all previous reconstructions by as much as 0.09 eV/(l x l)ce//. This indicates that the 
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adatom / restatom reconstruction is the most efficient way of saturating the highly unstable 
dangling bonds of the clean (111) surface. The alternative choice of a c(4 x 2) adatom / 
restatom geometry (using the two k-points obtained from the refolding in the c(4 x 2) ISBZ 
of the four ones used in the above (2 x 1) reconstruction calculations) confirms this result, 
yielding an energy 16 meV/(l x l)cell lower than the (2 x 2) choice. This energy difference 
is however comparable with our overall energy resolution. Similarly, our calculation does 
not rule out the existence of more complex surface reconstructions. It does however suggest 
that, if this is the case, then the adatom / rest atom mechanism, either in a (2 x 2) or in a 
c(4 X 2) arrangement, is very likely to constitute their basic building block, as is the case in 
Sicl and Gel. Moreover, the advantage over the vr-bonded chain reconstruction is larger in 
a-Sn, confirming the Si-Ge trend. 

In Table ^ we report the relaxed atomic coordinates of the (2 x 2) adatom-restatom 
reconstructed surface. Due to the symmetry constraint, the adatom and restatom only 
relax in the z direction. The in-plane distance between the rest atom and its three neighbors 
shrinks by about 10 % with respect to the ideal bulk value. The rest atom moves outward 
from its initial bulk-like position by 0.8 A , again a very large relaxation if compared with 
values in Si (0.3 A ), and in Ge (0.55 A ). Due to this relaxation, the rest atom has bond 
angles of 94.0°, a value closer to total s, p dehybridization (90°) than to the original sp^ 
one (109°), and bond lengths of 2.84 A (bulk bond length is 2.79 A ). In turn, the adatom 
bond lengths with the three neighboring first-layer atoms are 2.97 A , larger by 7 % than 
the bulk bond length, and the adatom bond angles are 93.4°. The second-layer atom under 
the adatom is pushed downward. Although its final distance from the adatom (2.99 A ) 
is comparable with the other adatom bond length, there is no bond-like accumulation of 
electronic charge between them, as can be seen in the total charge density reported in Fig. 
1^. Rather, most of the adatom charge is transferred to the restatom dangling bond. This 
charge transfer is not directly visible in the total charge density of Fig. |^, but shows up 
clearly in the band structure of Fig. ^. The half filled surface-localized band typical of 
unsaturated dangling bonds (see Fig. is here split into a lower filled band with prevailing 
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rest-atom character, and an upper empty band manly localized onto the adatom. This charge 
transfer is known to be the fundamental mechanism for the stabilization of the adatom / 
rest atom reconstruction in Si and Gei^l, and does appear to be so even in Sn. A glance at 
the charge density associated with the lower and upper surface bands (see Fig. ||) confirms 
this expectation. The charge associated with the lower band, shown in the panel (a) of Fig. 
H, has a strong rest-atom character and extends very little into the bulk, thus confirming its 
picture of saturated dangling bond. On the contrary, the upper unoccupied band (panel (b) 
of Fig. has a strong adatom character, but extends largely below the adatom, and might 
therefore be thought to as a band of so called "floating bonds" il associated with both the 
adatom and the five-fold coordinated atom immediately underneath. 

The optimal atomic positions and electronic band structure of the c(4 x 2) 
adatom/restatom reconstruction are given in Table ^ and Fig. 0. Differences in the 



atomic relaxations and charge densities are negligible with respect to the (2 x 2) case, and 
will not be re-discussed. 

We have also verified that the site for the adatom is not energetically favored, the 
energy of this surface being higher than the one with adatoms in T4 sites by about 70 
meV/(l X l)cell. Its optimal atomic coordinates are given in Table |V11. 



VII. METALLIC OVERLAYERS 

The presence of a metallic phase (/?-Sn) energetically very close to the a phase, suggests 
that an insulator-to-metal transition might take place at the surfaces of a-Sn, for example, 
through the formation of a thin /?-Sn metallic overlayer. Such was proposed to be the case, 
e.g., in gallium0, where based on calculations similar to these presented here a metallic 
bilayer of Ga-III was predicted to stabilize the surface of a-Ga better than any other recon- 
struction. In our case, however, we encounter a problem, since no low-index /5-Sn planes 
appear to match epitaxially the a-Sn (111) lattice. In particular, we have calculated the 
surface energy of several structures obtained by covering the a;-Sn(lll) surface with strained 
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epitaxial (100), (110), (111), and (221) planes of /3-Sn. None of these surfaces gave a sur- 
face energy comparable, let alone lower, than the ideal (111) surface, all our results being 
higher with respect to the (2 x 2) adatom / restatom reconstruction by more than 0.45 
eV/(l X l)cell. However, this negative result cannot of course be taken as a guarantee that 
metallic overlayers do not form at the a-Sn surface, since other more complex configura- 
tions, beyond our present factory, might have to be considered. The situation with respect 
to metallization remains therefore open. 

An indirect indication against metallization is however provided by surface energies. 
One can expect metallization to be favored, in fact, in cases where the metal has a lower 
surface free energy. In our case, however, the T=0 K energy of reconstructed a-Sn (111), 
Esurf ~ 540 mJ/m'^, is substantially lower than the free energy measured for /5-Sn0, namely 
670 mJ /m?. This suggests that metallization might not take place in the ground state, while 
it could be attained after deconstruction. 

VIII. DISCUSSION AND CONCLUSIONS 

From this ab-initio study, we predict that the (111) surface of a-Sn should be unstable 
in its unreconstructed form, and can be stabilized by various types of reconstructions. 

Among the reconstructions not involving adatoms, the (2x1) vr-bonded chain model is 
energetically best, and is predicted to display a gigantic buckling above 1 A in magnitude. 
This is in line with the trend towards increasing buckling in going from diamond to Si to 
Ge (111). 

In the presence of adatoms, the adatom-restatom pair reconstruction is found to be 
energetically best, and is thus the strongest candidate for the true ground state of this 
surface. The relative energy gain over the vr-bonded state is larger than in Si and Ge, again 
in line with trends. The actual optimal periodicity of the adatom-reconstructed a-Sn (111) 
is however difficult to predict (as is the case in Si and Ge) because of many subtle factors, 
including possible coexistence and competition between the two different basic units, namely 
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(2 X 2) and c(4 x 2). 

A search for a possible metallic state of a-Sn (111), suggested by the proximity of /?-Sn 
in the phase diagram, was so far inconclusive, and is disfavored by surface energy consider- 
ations. 

The beginning of actual experiments on a-Sn (111) is at this stage highly desirable, as we 
are not aware of any data so far. The present study makes rather strong predictions, both 
at the qualitative and the quantitative level. Qualitatively, we expect 7r-bonded chains to 
dominate on the adatom-free (e.g., cleaved) surface, and adatom-restatom reconstructions 
to appear after annealing, very much as for Si and Ge (111). At the quantitative level, 
we predict that buckling distortions and the associated dipole moments connected with 
intra-surface ionic charge transfers should be larger than in Si and Ge. This effect ought 
to be relevant to, e.g., surface vibrational spectroscopy. The corresponding band splittings 
should finally be observable in electron spectroscopies, including photoemission and optical 
absorption, particularly to optical conductivity ior q ^ 0. 

It is hoped that these results will stimulate newer efi^orts towards an understanding of 
the CK-Sn (111) surface, where hardly any data can presently be found. 
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TABLES 
TABLE I. 





C 


Si 


Ge 


2x1 TT-bonded chains 


dimerized (*) (?) 


buckled 


more buckled 


adatom/ rest atom 


not found 


7x7 DAS model (*) 


c(2 X 8) adatom model (*) 


2x1 buckled 


not found 


saddle-point 


saddle-point 


metallic 


graphitic T > 2700K (?) 


T > 1500K (?) 


T > 1050K 



(*) Stable structure at T = 0. 
(?) Pending confirmation. 



TABLE IL Calculated surface energies, absolute and relative, of different optimized reconstruc- 
tions for a-Sn (111) surface. 



Structure 


-^surf 


AE 


vert, relax, or buckling^ 




[eV/(lxl cell)] 


[mJ/m?] 


[eV/(lxl cell)] 


(A) 


ideal 


0.940 


837 


0.000 





fully relaxed 


0.933 


831 


-0.007 


0.02 


2x1 buckled (Haneman) 


0.732 


652 


-0.217 


1.23 


2x1 TT-bonded chain 


0.697 


621 


-0.243 


1.15 


2x2-adatom {H^) 


0.696 


620 


-0.244 


0.54 (restatom) 


2x2-adatom {T^) 


0.626 


558 


-0.314 


0.81 (restatom) 


c(4x2)-adatom (T4) 


0.610 


543 


-0.330 


0.74 (restatom) 


metallized 


> 1.06 


944 


> 0.121 





^Vertical relaxation is relative to initial bulk-like positions, and buckling is difi^erence between 
relaxation of two atoms. 
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TABLE III. Ideal, and optimized atomic positions of the a-Sn (111) (2 x 1) buckled (Haneman) 
surface. In the rectangular supercell, coordinates are given by r = Ciai+C2a2+C3a3, where is de- 
fined in the conventional cubic coordinate system as ai = (ao/2)(— 1,2, —1), a2 = (ao/2)(— 1,0, 1), 
^3 = ao(l) 1) 1) and gq (=6.446 A ) is the lattice parameter. 



Atom 




ideal 






Optimal 




no. 


Cl 






Cl 


C2 


C3 


1 


.[}[}[} 


.(JUU 


.uuu 


.UU9 


.UUU 


-.U45 


o 
Z 


.5111) 


.5(JU 


.UUU 


.4oo 


.oUU 


.Udz 


6 


.167 


.500 


.083 


.132 


.500 


ATT 

.077 


4 


.667 


.000 


.083 


.692 


.000 


.077 


5 


.167 


.500 


.333 


.169 


.500 


.330 


6 


.667 


.000 


.333 


.667 


.000 


.332 


7 


.333 


.000 


.417 


.334 


.000 


.415 


8 


.833 


.500 


.417 


.836 


.500 


All 


9 


.333 


.000 


.667 


.334 


.000 


.665 


10 


.833 


.500 


.667 


.833 


.500 


.667 
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TABLE IV. Initial (Pandey's), and optimized atomic positions of the a-Sn (111) (2 x 1) 
vr-bonded surface. In the rectangular supercell, coordinates are given by r = ciai + C2a2 + c^a.^, 
where aj is defined in the conventional cubic coordinate system as ai = (ao/2)(— 1, 2, — 1), 
^2 = (ao/2)(— 1,0, 1), a3 = ao(l, 1, 1) and oq (= 6.446 A ) is the lattice parameter. 



Atom 
no. 


Cl 


Pandey's*^ 


C3 


Cl 


Optimal 


C3 


-1 
1 


1 /I Q 

-.145 


.UUU 




-.lb / 


.000 


.UoU 


2 


-.019 


.500 


.023 


-.020 


.500 


-.043 




.610 


.500 


.106 


.278 


.500 


.100 


4 


.519 


.000 


.106 


.484 


.000 


.077 


5 


.167 


.500 


.333 


.178 


.500 


.341 


6 


.667 


.000 


.333 


.657 


.000 


.317 


7 


.333 


.000 


.417 


.336 


.000 


.431 


8 


.833 


.500 


.417 


.831 


.500 


.400 


9 


.333 


.000 


.667 


.337 


.000 


.678 


10 


.833 


.500 


.667 


.834 


.500 


.654 



'^For corresponding bulk-like positions, see the ideal case of Table |TT|. 
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TABLE V. Initial, and optimized atomic positions of the a-Sn (111) (2x2) adatom / re- 
statom (T4 configuration) reconstructed surface. In the hexagonal supercell, coordinates are given 
by r = ciai + 02^2 + c^a^, where is defined in the conventional cubic coordinate system as 
ai = (oo/2)(-l,2, -1), a2 = (ao/2)(-l,0, 1), as = (ao/3)(l, 1, 1) and oq (= 6.446 A ) is the lattice 

paramolcr. 



Atom 
no. 


Cl 


Initial 

C2 


C3 


Cl 


Optimal 

C2 


C3 


1 


.000 


.000 


2.375 


.000 


.000 


2.509 


2 


-.167 


.500 


2.125 


-.158 


.475 


2.077 


3 


.333 


.000 


2.125 


.317 


.000 


2.077 


4 


-.167 


-.500 


2.125 


-.158 


-.475 


2.077 


5 


.333 


-1.000 


2.125 


.333 


-1.000 


2.343 


6 


.000 


.000 


1.875 


.000 


.000 


1.705 


7 


.500 


.500 


1.875 


.485 


.544 


1.935 


8 


.000 


-1.000 


1.875 


.030 


-1.000 


1.935 


9 


.500 


-.500 


1.875 


.485 


-.544 


1.935 


10 


.000 


.000 


1.125 


.000 


.000 


.963 


11 


.500 


.500 


1.125 


.500 


.501 


1.175 


12 


.000 


-1.000 


1.125 


.000 


-1.000 


1.175 


13 


.500 


-.500 


1.125 


.500 


-.501 


1.175 


14 


.167 


.500 


.875 


.175 


.525 


.864 


15 


.667 


.000 


.875 


.667 


.000 


.926 


16 


.167 


-.500 


.875 


.175 


-.525 


.864 


17 


.667 


-1.000 


.875 


.650 


-1.000 


.864 


18 


.167 


.500 


.125 


.170 


.510 


.110 


19 


.667 


.000 


.125 


.667 


.000 


.175 


20 


.167 


-.500 


.125 


.170 


-.510 


.110 



20 



21 .667 -1.000 .125 .660 -1.000 .110 
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TABLE VI. Initial, and optimized atomic positions of the a-Sn (111) c(4 x 2) adatom / re- 
statom (T4 configuration) reconstructed surface. In the rectangular supercell, coordinates are 
given by r = ciai + C2a2 + c^a^, where is defined in the conventional cubic coordinate system as 
ai = (ao/2)(-l,2, -1), a2 = (ao/2)(-l,0, 1), as = (ao/3)(l, 1, 1) and ao (=6.446A ) is the lattice 

paramolcr. 



Atom 
no. 


Cl 


Initial 

C2 


C3 


Cl 


Optimal 

C2 


C3 


1 


.000 


.000 


2.375 


.004 


.000 


2.527 


2 


.167 


.500 


2.125 


.158 


.469 


2.083 


3 


-.333 


.000 


2.125 


-.313 


.000 


2.102 


4 


.167 


-.500 


2.125 


.158 


-.469 


2.083 


5 


-.333 


-1.000 


2.125 


-.341 


-1.000 


2.325 


6 


.000 


.000 


1.875 


-.006 


.000 


1.716 


7 


-.500 


.500 


1.875 


-.489 


.525 


1.941 


8 


.000 


-1.000 


1.875 


-.027 


-1.000 


1.945 


9 


-.500 


-.500 


1.875 


-.489 


-.525 


1.941 


10 


.000 


.000 


1.125 


.001 


.000 


.975 


11 


-.500 


.500 


1.125 


-.501 


.501 


1.178 


12 


.000 


-1.000 


1.125 


.005 


-1.000 


1.182 


13 


-.500 


-.500 


1.125 


-.501 


-.501 


1.178 


14 


-.167 


.500 


.875 


-.174 


.525 


.880 


15 


-.667 


.000 


.875 


-.650 


.000 


.874 


16 


-.167 


-.500 


.875 


-.174 


-.525 


.880 


17 


-.667 


-1.000 


.875 


-.665 


-1.000 


.907 


18 


-.167 


.500 


.125 


-.169 


.509 


.127 


19 


-.667 


.000 


.125 


-.661 


.000 


.123 


20 


-.167 


-.500 


.125 


-.169 


-.509 


.127 
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21 -.667 -1.000 .125 -.668 -1.000 .149 
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TABLE VII. Initial, and optimized atomic positions of the a-Sn(lll) (2 x 2) adatom / re- 
statom {H^ configuration) reconstructed surface. In the hexagonal supercell, coordinates are given 
by r = ciai + C2a2 + csSls, where aj is defined in the conventional cubic coordinate system as 
ai = (ao/2)(-l,2,-l), a2 = (ao/2)(-l, 0, 1), as = (ao/3)(l, 1, 1) and oq (= 6.446 A ) is the lattice 

paramolcr. 



Atom 
no. 


Cl 


Initial 

C2 


C3 


Cl 


Optimal 

C2 


C3 


1 


.167 


.500 


2.375 


.167 


.500 


2.525 


2 


-.167 


.500 


2.125 


-.142 


.500 


2.090 


3 


.333 


.000 


2.125 


.321 


.037 


2.090 


4 


-.167 


-.500 


2.125 


-.167 


-.500 


2.271 


5 


.333 


-1.000 


2.125 


.321 


-1.037 


2.090 


6 


.000 


.000 


1.875 


-.015 


-.046 


1.848 


7 


.500 


.500 


1.875 


.531 


.500 


1.848 


8 


.000 


-1.000 


1.875 


-.015 


-.954 


1.848 


9 


.500 


-.500 


1.875 


.500 


-.500 


1.984 


10 


.000 


.000 


1.125 


.000 


-.001 


1.102 


11 


.500 


.500 


1.125 


.501 


.500 


1.102 


12 


.000 


-1.000 


1.125 


.000 


-.999 


1.102 


13 


.500 


-.500 


1.125 


.500 


-.500 


1.199 


14 


.167 


.500 


.875 


.167 


.500 


.861 


15 


.667 


.000 


.875 


.663 


-.012 


.883 


16 


.167 


-.500 


.875 


.175 


-.500 


.883 


17 


.667 


-1.000 


.875 


.663 


-.988 


.883 


18 


.167 


.500 


.125 


.167 


.500 


.105 


19 


.667 


.000 


.125 


.665 


-.006 


.133 


20 


.167 


-.500 


.125 


.171 


-.500 


.133 
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21 .667 -1.000 .125 .665 -.994 .133 
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FIGURES 



FIG. 1. Electronic structure of bulk a-Sn. The zero in energy corresponds to the Fermi level. 
FIG. 2. Brillouin zone and electronic structure of bulk /?-Sn. The zero in energy corresponds 
to the Fermi level. 

FIG. 3. Surface electronic structure of the ideal Q;-Sn(lll) surface reported along 
high-symmetry lines of the (1 x 1) hexagonal irreducible Brillouin zone. Shaded areas correspond 
to surface-projected bulk states, while thicker lines correspond to surface states. The surface Bril- 
louin zone is given in the inset. Note that the dangling bond surface state crosses the fermi level 
Ep. 

FIG. 4. Electron density contour of the highest occupied state of the ideal (1 x 1) surface at 
K, on the (110) plane passing through top atoms. Full circles correspond to Sn atoms, and thicker 
straight lines to bonds among Sn atoms. Contour lines are separated by 0.0005 (a.u.). Note the 
dangling bond character of this state. 

FIG. 5. Electron density contours of the highest occupied state ( panel (a) ) and of the lowest 
unoccupied state ( panel (b) ) in the (2 x 1) vr-bonded chain reconstruction at J, on the (110) plane 
passing through the up (down) atom. In panel (b) atoms are labelled according to the optimal 



positions of Table IV. Contour lines are separated by 0.0005 (a.u.). 

FIG. 6. Surface electronic structure of the (2 x 1) vr-bonded chain reconstructed surface re- 
ported in the (2 x 1) rectangular irreducible Brillouin zone (see inset). Shaded areas correspond to 
surface-projected bulk states, while thicker lines correspond to surface states. The original dangling 
bond is now split into two. 

FIG. 7. Geometry, and electron density contour of the (2x2) adatom / restatom reconstruction, 
on the plane passing through the adatom and the rest atom. Contour lines are separated by 0.005 
(a.u.). Note the strong outward relaxation of the restatom, and the strong inward relaxation of 
the second and third layer atoms beneath the adatom. (R: restatom. A: adatom) 
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FIG. 8. Surface electronic structure of the (2 x 2) adatom / restatom reconstructed surface 
reported in the (2 x 2) hexagonal irreducible Brillouin zone (see inset). Shaded areas correspond 
to surface-projected bulk states, thicker lines correspond to surface states, dotted lines to surface 
resonances. The ~ 1 eV splitting of the surface state reflects the adatom-restatom electron transfer. 

FIG. 9. Electron density contours of the highest occupied state ( panel (a) ) and of the lowest 
unoccupied state ( panel (b) ) in the (2 x 2) adatom / restatom reconstruction at K, on the same 
plane as in Fig. |^ In panel (b) atoms are labelled according to the optimal positions of Table 0. 
Contour lines are separated by 0.0004 (a.u.). Note the strong restatom / adatom characters, with 
a larger penetration of the adatom empty state. 

FIG. 10. Surface electronic structure of the c(4 x 2) adatom / restatom reconstructed surface 
reported in the c(4 x 2) rectangular irreducible Brillouin zone (see inset). Shaded areas correspond 
to surface-projected bulk states, thicker lines correspond to surface states, dotted lines to surface 
resonances. 
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